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Summary 

This Quarterly Report addresses the following topics related to the dynam- 
ics of the 4-mass tethered system: 

The development of damping algorithms for damping the out-of -plane libra- 
tion of the system and the interaction of the out-of -plane control with the 
other degrees of freedom. 

The development of environmental models to be added to the dynamics simula- 
tion computer code. The environmental models are specifically a new drag rou- 
tine based on the Jacchia’s 1977 model, a J 2 model and an accurate thermal model 
of the wire. 
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Figure Captions 


Figure 1. 


Figure 2 a~b. 


Figure 3a-b. 


Figure 4a-i. 


Figure 5a-b. 


Figure 6a-r. 


Figure 7a-b. 


Figure 8. 


Figure 9a-d. 


Figure 10. 


Figure 11. 


In-plane and out-of-plane libration angles vs. time for a case 
in which the out-of -plane damping system is disactivated . The 
gains of the control laws are: K^=2 , K^=0 . 

Dynamic response vs. time for a case in which a modified 
(Jp = constant) yo-yo control law is adopted for out-of- 
plane libration control while the in-plane control system 
is disactivated. The gains of the control laws are: 

K*=0, K^=4. 

The same as in Figure 2 except that this time the in -plane 
control system is activated. The gains of the control laws 
are: K^=l , K^=4 . 

The same as in Figure 2 (the in-plane control system is 
disactivated) . The gains of the control laws are: 

K*=0 , K*>=10. 

The same as in Figure 3. The gains of the control laws 
are: K e =0.1, K^=10. 

Dynamic response vs. time for a case in which the 
standard yo-yo control law (variable <p) is adopted 
for out-of-plane libration control . The in-plane control 
system is also activated. The gains of the control laws 
are: K 0 =O . 2 , K*=10. 

Front-view and side-view respectively of the system under 
the conditions of Figure 6. The snapshots are every 100 
sec for 36,000 sec. The x-axis scale is largely expanded. 

Atmospheric density vs. altitude for different exospheric 
temperatures. The curves are obtained by means of an 
analytical fit of the Jacchia's 1977 model. 

Atmospheric density vs. time and exospheric temperature 
vs. time over one orbit at 450 km altitude and 28°.5 
inclination. Figure 9a-b show the results for the sun at 
the Spring Equinox and Figure c-d for the sun at the Summer 
Solstice . 

Acceleration components in the inertial reference frame 
due to the J 2 gravity term vs . time over one orbit 
at 450 km altitude and 28°.5 inclination. 

Acceleration components in the orbiting reference frame 
due to the J 2 gravity term vs. time for the same 
conditions of Figure 10. 
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1.0 INTRODUCTION 

This is the sixth Quarterly Report submitted by SAO under contract NAS8- 
36606, "Analytical Investigation of the Dynamics of Tethered Constellations in 
Earth Orbit (Phase II)," Dr. Enrico C. Lorenzini, PI. This report covers the 
period from 22 June 1986 through 21 September 1986. 


2.0 TECHNICAL ACTIVITY DURING REPORTING PERIOD AND PROGRAM STATUS 

2.1 Algorithms For Damping Out-Of-Plane Librations 

Damping of out-of-plane librations is achieved by controlling the lengths 
of the three tether segments connecting the four platforms. The out-of-plane 
librational dynamics of a multi-mass system can be modeled, with good approxima- 
tion, by reducing the system to a two-mass system if the bowing of the former is 
small. The simplified equation of motion that describes the unperturbed, small 
out-of-plane librations of a two-mass tethered system is as follows: 

l 2 <p + 2 ti<f> + 4£ 2 nv = 0 (2.1.1) 


Eor a two-mass tethered system the energy dissipated per cycle per unit mass is 
therefore given by ( i ~ t Q for St « t 0 ) ' 



where r ^ is the period of the out-of-plane libration equal to tt/O. From simple 
physical considerations and also from equation (2.1.1') it is immediately evi- 
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dent that the energy dissipated during an out-of-plane damping cycle does not 
depend upon the orbital rate 0. The out-of-plane libration is in fact very 
weakly coupled with the longitudinal degrees of freedom, consequently tension 
control laws, which rely on such coupling, are poorly effective in damping out- 
of-plane librations. A simple mechanical equivalent of the out-of-plane dynam- 
ics of a tethered system is a pendulum of variable length on the ground (the 
pendulum has negligible Coriolis forces due to the Earth rotational rate and the 
librational frequency is approximately constant if the length variations are 
small) . The yo-yo technique developed in reference [1] and [2] are perfectly 
suitable for damping the out-of-plane librations of a tethered orbital system. 
The same technique is on the contrary surpassed by the S-type [3] control law 
for damping the in-plane libration because the latter control law exploits the 
Coriolis forces generated by the orbital rate Q. In order to clarify this point 
we can use equation (2.1.1*) for computing the energy dissipated by substituting 
into that equation various control laws. In the case of small damping we can 
assume that the libration over one cycle is approximately sinusoidal. We have 
therefore 


<p ~ y?sin (2flt) (2.1.2) 

<p ~ 2C2y?cos (2ftt) (2.1.3) 

An ideal yo-yo control law (bang-bang type) can be expressed in formulae as 
follows : 


l/e o = (1 +K v <p) if sign (<£>) = sign ( <p ) 

l/l 0 = (l-K^p) if sign (<p) ^ sign ( <p ) 


(2.1.4) 



Page 8 


where <p is the amplitude of the out-of-plane libration. We can rewrite equa- 
tion (2.1.1*) as follows: 





v> 

d <p 



d<p 


(2.1.5) 


where is the out -of -plane component of the damping torque per unit mass and 
L ^ is the out-of-plane component of the angular momentum. By making use of 
equation (2.1.4) we finally have: 


— — 2 ( l a** lain ) IPmax — “ 320 3 ^p 3 K^ £ 0 


( 2 . 1 . 6 ) 


The fractional damping of the <p degree of freedom can be obtained by equating 
the kinetic energy variation for that degree of freedom to the energy dissipated 
per cycle (per unit mass) given by equation (2.1.6). We have: 


ATp — Ed , tp 


(2.1.7) 


or 


2 tin 2 [(£-A£) J -£ 2 ] 32 n J ^ 3 K^* 


( 2 . 1 . 8 ) 


After defining: 


St — l ■ ax ^min — 2Kip<pt 0 


(2.1.9) 


we obtain : 


A <p . St 

<P t o 


(ideal yo— yo, out - of - plane) 


( 2 . 1 . 10 ) 
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An ideal yo-yo control law, as expressed by equation (2.1.4) is unfeasible be- 
cause it requires an infinite tether velocity. In actuality the ideal yo-yo 
control law is approximated by a sinusoidal control law (we call it smooth yo- 
yo) that can be defined as: 


t/t 0 = [1 + K^sin (40t)] (2.1.11) 

By substituting equation (2.1.11), its derivative and equation (2.1.3) into 
equation (2.1.1') we obtain: 

( 2 . 1 . 12 ) 

After making use of equation (2.1.7) and of equation (2.1.9) we finally obtain: 
AS 8 1 

~ 7T — (smooth yo-yo, out - of —plane) (2.1.13) 

<p to 

By using the smooth yo-yo instead of the ideal yo-yo we have therefore a de- 
crease in damping effectiveness of roughly 20%. 

A modification of the smooth yo-yo control law can be obtained by assuming 
that ip in equation (2.1.11) is constant and equal to the amplitude of the first 
out-of-plane libration cycle Xp\. Consequently equation (2.1.13) is modified to 

Si ir (2.1.13') 

V V t a 

The fractional damping is therefore increasing with the decrease of ip. This 
result seems theoretically remarkable. In practice, however, the control of 
out-of-plane librations excites in-plane librations because the latter are 
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strongly coupled with tether length variations. This spill-over effect pumps 
the in -plane libration amplitude. By adopting the modified smooth yo-yo control 
law the spill-over does not decrease with <p\ the in-plane libration builds up. 
The in-plane control consequently tries to damp out the in-plane libration and 
by doing so impairs the effectiveness of the out-of -plane libration damping. In 
actuality therefore the result shown by equation (2. 1.13*) is not achievable. 

We want now to compare the in-plane libration damping with the out-of-plane 
libration damping. The simplified equation of motion that describes the unper- 
turbed, small in-plane librations of a two-mass tethered system is as follows: 

i 2 e + 2 ii(»-n) + 3 t 2 n 2 e = o (2.1.14) 

The energy dissipated over a cycle per unit mass is therefore given by ( t ~ i 0 
for St « t 0 ) : 


. r T B 

ii{ 6 -n) 6 dt a - 2 1 0 i( 0 -a)ddt. ( 2 . 1 . 15 ) 

J o 

where tq is the in-plane libration period equal to 2ir/ ( \/3~ft) . The equation, 
similar to equation (2.1.7), applicable to in-plane librations can be obtained 
by multiplying equation (2.1.14) by 0 and intergrating from o to Since 

d 0 2 /dt = 209 and t ~ t 0 we obtain: 

/ Tff [T0 mT0 

d 6 2 + 2t 0 j i( 6 -n) 6 dt + 3 n 2 £o J OdO = o (2.1.16) 

Equation (2.1.16) can therefore be written as follows: 



AT* r#1 = Ed,* 


(2.1.17) 
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where AT$ stands for the kinetic energy of the 0 degree of freedom relative to 
the LH-LV reference frame. 

The smooth yo-yo control law for damping the in-plane libration is ex- 
pressed as follows: 


t/t„ = [l+K 9 0sin(2 v^Tnt)] 


(2.1.18) 


By following the same procedure adopted for computing the out-of-plane frac- 
tional damping but using equation (2.1.17) instead of (2.1.7) we arrive at the 
following result: 


— ~ 7C — (smooth yo-yo, in -plane) (2.1.19) 

$ to 

This result stems from the fact that by using a yo-yo type control law we do not 
exploit the Coriolis forces generated by the orbital rate 0. The fractional 
damping turns therefore to be the same for in-plane librations and out-of-plane 
librations; it would be completely different had we used a tension control law. 
The conclusion just reached prompted us to design the S-type control law for 
damping the in-plane libration. The S-type control law is formulated as fol- 
lows : 


i/to = (1-K*0) (2.1.20) 


The associated fractional damping, in the case of light damping, is given by: 


Afl ^ y St 

T ~ \/ r 3~ T 0 e 


( S — type , ir. — plane) 


( 2 . 1 . 21 ) 
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The most important feature is that the S-type control law provides a fractional 
damping that is independent of the amplitude 1 . Heavy damping can be easily 
achieved with this type of control law even at small amplitude 1. In the case 
of heavy damping, however, equation (2.1.21) provides only a first order approx- 
imation of the fractional damping because the in-plane libration is no longer a 
pure sinusoid over a cycle; equation (2.1.21) gives an overestimation of the 
fractional damping in that case. 

We have also computed the fractional damping for out-of-plane librations 
achieved by using other non linear control laws. Table 1 summarizes the frac- 
tional damping achievable vs. the various control laws shown in column one for 
the out-of-plane dynamics. Table 2 shows similar results for the in-plane dy- 
namics . 


Table 1 

Out-of-Plane Libration Damping by Means of Tether-Length Control Law 

Contro l Lau Mathematical Formulation Fractional Damping 

Ideal Yo-Yo t/t 0 ={ 1+K„^) if sign (<p) = (<p) ^ ~ 4 — 

to 

t/t 0 =(l-Krp) if sign (<p) ^ sign (<p) 


Smooth Yo-Yo 

l/l 0 = [l+K„^sin(4nt)] 

A <p 6 £ 

~ n — 



<P 



Modified 

t/t 0 = [1+Kp^iSin (4flt) ] 

A <p £>i 

61 

Smooth 


<P 

<P 

to 


Yo-Yo 


A <p 
<P 


IT 


61 

to 


<p<p Control 
Law 


t/lo = [1+K r<pp] 
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Table 2 


In -Plane Libration Damping by Means of 

Tether-Length Control Laws 

Control Law 

Mathematical Formulation 

Fractional Damping 

Ideal Yo-Yo 

l/l 0 -( 1+K*0) if sign (0) = (0) 

t/t 0 = (1-M) if sign (0)^(0) 

a! . st 

/N/ A 

1 ~ t o 

Smooth Yo-Yo 

l/t 0 = [1+K$0sin (2 y/$~ Ot) ] 

A0 St 

— ^ 7T — 

$ t 0 

S -Type 

t/t 0 = (1-K e 0) 

A0 ~ ir St 
o ~ \/T~ tj 

References to 

Section 2.1 



1. Baker W.P. et al . "Tethered Subsatellite Study," NASA 
TM X-73314 , March, 1976. 

2. Spenser T.M. "Atmospheric Perturbation and Control of 

a Shuttle/Tethered Satellite," Proceedings of the 8 th IEAC 
Symposium, Automatic Control in Space, 1980. 
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2.2 Simulation Runs With/Without Out-Of -Plane Libration Damping Algorithms 

The previous section dealt with the yo-yo damping technique for a two-mass 
tethered system. In this section we generalize it to a four-mass tethered sys- 
tem. As already shown for the in-plane control system in Quarterly Report No. 5 
Section 2.4 the generalization is readily implemented by controlling the three 
tether segments proportionally to their unstretched lengths. In formulae 

^lc = £oi (1 — K $6 +K € i 1 €ii/£oi + ^dlol) 

^ 2c = Zo 2 (1 — K$B + oa — K € ii€ii/£o2 + ^dlol) (2.1.22) 

^ 3c = ^03 (1 — K$& — 02 // 0 3 + ^DLOl) 

Equations (2.1.22) are similar to equations (2.4.3) of Quarterly Report No. 5 
except for the term £ D lol that is the tether length variation per unit length for 
the out-of-plane libration control. The variable £dlol depends on the particular 
control law adopted. 

Several simulation runs have been done in order to derive a reasonable 
technique for damping the out-of-plane oscillations. The smooth yo-yo control 
law has been extensively tested either in the standard (variable p) or in the 
modified version (constant p) . The in-plane libration angle B and the out-of- 
plane libration angle p are shown in Figure 1 for a case in which the in-plane 
libration control is activated (K#=2) while the out-of-plane libration control 
is off (K^=0) . Notice that the in -plane libration is very well damped and that 
the in-plane control does not spill over the out-of-plane libration. The vice- 
versa is not true. The second simulation was run by adopting K^=0 (in-plane 
control off) and in the modified smooth yo-yo control law. Figure 2a 

shows the out-of-plane angle vs. time. The fractional damping (for light damp- 
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Ing: fractional damping ~ logarithmic decrement = A; damping ratio ~ — ) of 

2ir 

the angle <p is approximately 11%. Figure 2b shows the in-plane angle and the 
out-of-plane angle vs. time for the same case. This figure clearly shows the 
spill of the out-of-plane control over the in-plane libration. The amplitude of 
the in-plane angle is building up. A third simulation was run in order to 
address the issue of the back action of the in-plane control over the out-of- 
plane damping. The control law adopted for the out-of-plane libration is the 
same of the previous run with the same gain k^ s= 4. The in-plane control is 
activated with a gain ko = 1 in the S-type control law. Figure 3a shows the 
out-of-plane angle vs. time for a time duration of 5000 sec. Because of the 
activation of the in-plane control the out-of-plane damping is reduced from 11% 
to 6.5%. The in-plane amplitude growth due to the spill-over is, on the con- 
trary, reduced from 2° to 1° as shown in Figure 3b. In the fourth simulation 
run K#=0 and K^—10 . This simulation has been run for 10,000 sec to show better 
the build-up of in-plane librations. Figure 4a depicts the out-of-plane libra- 
tion angle vs. time. A fractional damping ~ 26% has been achieved. The out- 
of-plane libration phase-plane, depicted in Figure 4b, stresses this result. 
Figure 4c shows the in -plane libration angle build-up (together with the out-of- 
plane angle vs. time) while Figure 4d depicts the phase-plane of the in-plane 
libration. The fractional tether length variation 6i/t 0 for out-of-plane libra- 
tion control is shown in Figure 4e (in the case of the modified smooth yo-yo, 
6t/t 0 has a constant amplitude). By substituting the values of 6t/l 0 = 0.12 
(remember that ^ is a peak-to-peak variation) into equation (2.1.13*) for the 
first cycle (when <pi/<p = 1) we obtain a fractional damping A 6/6 = 38% that is 
greater than the measured one (26%) . The difference can be attributed to two 
different effects: 1) Equation (2.1.13*) and all the other similar equations 

are valid for light damping and they overestimate the damping if otherwise. 2) 
The tether elasticity deforms the eight-shaped damping cycle invoked by equation 
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(2.1.11). This latter effect is readily noticeable in Figure 4f that shows the 
trajectory of mass no. 1 with respect to the orbiting reference frame (located 
approximately at the system C.M.). For completeness Figure 4g and 4h show the 
in -plane and out-of -plane components respectively of the distances of the mid- 
dle-masses from the straight line through the end-masses (lateral deflections) . 
Figure 4i depicts the moduli of such distances. In the fifth simulation the 
same gain k^ = 10 of the previous simulation was adopted for the out-of -plane 
libration control while K$ was set equal to 0.1 in an attempt of limiting the 
in-plane oscillation growth without impairing the damping of out-of-plane libra- 
tions . Figure 5a shows the out-of-plane libration angle vs. time. The frac- 
tional damping reduction is almost negligible: from 26% to 25%. Figure 5b 
shows the in-plane libration angle vs. time and the out-of-plane angle vs. time. 
The in-plane amplitude growth is reduced compared to the previous run. All the 
other features of this run are similar to the previous one. Because of the 
results obtained by adopting the modified yo-yo control law we decided to switch 
to the standard yo-yo (variable <p) for the following simulations in order to 
reduce the spill-over due to the out-of-plane libration control. 

The sixth simulation run adopts 1^=10 in the standard yo-yo control law 
while Ke is set equal to 0.2. This simulation has been run for 36,000 sec. 
Figure 6a shows the angle <p vs. time. The fractional damping is, at the begin- 
ning, around 17%. The fractional damping starts decreasing past 15000 sec. The 
variation is, however, non-monotonic: the fractional damping goes almost to 
zero around 30000 sec and goes up again to 6% at 32000 sec. Various sources can 
contribute to the deviation of the fractional damping behavior from equation 
(2.1.13): 1) the tether length variation for yo-yo control, as expressed in 
equation (2.1.11) , has no frequency feedback; the controller goes slowly out of 
phase with the <p oscillation because of the frequency variation due to damping. 
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This phenomenon, however, is not very important in the present simulation as it 
clearly appears from Figure 6b that shows the tether length variation for out- 
of -plane libration damping (£dlol) vs. the out- of -plane angle. The phase differ- 
ence is fairly constant as it can be inferred by the crossing points of the 
eight-shaped curves close to the 0 vertical line. 2) The system is probably 
settling down to a state in which energy, through exchange among the various 
degrees of freedom, is flowing into the least damped degrees of freedom. This 
conclusion can be preliminarily inferred from Figure 6c that shows the angular 
foot-print of the system (out-of -plane libration vs. in-plane libration). At 
approximately the same time when the damping starts decreasing the system gets 
trapped in an eight-shaped angular foot-print. 

Figure 6d shows the in-plane angle vs. time along with the out-of-plane 
angle vs. time. The in-plane angle is clearly forced by the tether length 
variation for the control of out-of-plane librations. By adopting the standard 
smooth yo-yo control law, however, the spill-over decreases with time. Figures 
6e and 6f are the phase-planes for the <p and 0 oscillations respectively. Fig- 
ures 6g and 6h show the tether length variation per unit tether length for the 
out-of -plane and in-plane oscillation control respectively. Figure 6i and 6j 
show the in-plane and out-of-plane components respectively of the lateral 
deflections of the Space Station (solid line) and the G-laboratory (dotted 
line) . The in-plane component of the lateral deflection of the G-laboratory is 
excited at the beginning of the simulation by the out-of-plane control spill- 
over but it decreases as the spill-over decreases. The out-of-plane components 
of lateral deflections of both the Space Station and the G-laboratory increase 
slowly with time. Since there is no control over such components in all these 
simulation runs the energy is probably accumulating in these degrees of freedom. 
Figures 6k, Si and 6m depict the tether length variations for tether no. 1, no. 
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2 and no. 3 respectively. In order to give an idea of the acceleration fluctua- 
tions during the out-of-plane damping process Figure 6 t depicts the three compo- 
nents, in the orbiting reference frame, of the acceleration on board the Space 
Station. The residual acceleration is due to the fact that all the initial 
oscillations have not been damped out completely. The in-plane libration how- 
ever can be damped out very effectively by disactivating the out-of-plane yo-yo 
control and by increasing the gain of the in-plane control system. The follow- 
ing figures, therefore, should not be interpreted as the limit performance of 
the system. Figures 6p, 6q and 6r show the flight direction (ACH) , transverse 
(ACL) and vertical (ACH) component of the acceleration on board the G-laboratory 
vs. time. Because of the relatively large residual libration amplitudes the 
residual acceleration fluctuations on board the G-laboratory are larger than the 
acceleration fluctuations on board the Space Station. Finally Figures 7a and 7b 
show the expanded-scale (the expansion coefficient is different for the two 
figures) front-view and side-view respectively of the system with snapshots ev- 
ery 10 0 sec up to 36,000 sec. The residual oscillations of the system are 
clearly evident in these two figures. 


2.3 Environmental Models 

We have developed an entirely new package of subroutines which model the 
environment at the altitudes of interest to the Space Station. 

We have modeled the major perturbations acting upon a tethered system in 
the range of altitudes from 400 km to 600 km. The major perturbations are: 
drag, J 2 gravity term and thermal as dealt with in detail in the next subsec- 


tions . 
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2.3.1 Drag Model - 

As a density model we have adopted the version by C.E. Roberts Jr., refer- 
ence [1] , who developed an analytical model of upper atmosphere which provides 
atmospheric density as a function of exospheric temperature and altitude. The 
model is based on Jacchia’s 1970, taking also into account the '77 corrections 
[2] . The temperature equation has been modified in order to have a simple 
analytical expression of the temperature as a function of the height h. 

The atmospheric density is given by: 

6 


p(h) = Y, Midi(h) (2.3.1) 


where : 


1=1 


i 


constituents = N 2 , Ar , He, 0 2 , 0, H 


Mi 


molecular mass (gr/mole) 


di 

di (h) 

= 

number of particles of each constituent per cm 3 /Avogadro 1 s 
number 

di (h 125 ) 

= 

computed data at h = 125 km 


Tex 

= 

exospheric temperature 


<*i 


thermal diffusion coefficients 


T* 

= 

f (exospheric temperature) = A x + A 2 T EX + A 2 exp (A 4 T EX ) 

T (h) 

= 


(2.3.3) 

T 0 

= 

188 °K 


h x 

= 

125 km 


h 0 

= 

90 km 


l 

— 

to be chosen in order to produce a minimum in the absolute 
value of the percentage difference with the Jacchia's 
model for altitudes between 125 and 1000 km and for T E x 
between 600°K and 2000°K 
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R* 

7i 


g Q 

R 


effective earth's radius = 6356.766 km 
Mjg 0 R m / Tex ~~ T x \ / h x — hp \ 

R^T ex V T EX-ToA R »+ h x/ 

gravity acceleration on the ground = 9.80665 — 

s 2 

universal gas constant = 8.31432 J/°K-mole 


(2.3.4) 


The number of particles of atomic hydrogen d H is equal to zero for h < 50 0 km. 
As a result of the above mentioned computations Figure 8 shows the atmospheric 
density vs. altitude, in the range of interest to our system, for different 
exospheric temperatures. 


In order to obtain a more accurate density model we have improved the 
"static model" by taking into account the phenomena which affect the values of 
the global exospheric temperature [3] . We have modeled the effects of solar 
activity and diurnal variations but we have neglected the geomagnetic activity. 
Following these assumptions we can express the local exospheric temperature T L 
as follows: 


where : 


T l = T c (1 +R T sin m 0) ( 1 +R T 


cos m y? — sin n fl 
1 +R T sin m 0 


-§) 


(2.3.5) 


T c 


R x 


m 

n 


V 


<t> 


night minimum of the global exospheric temperature, 
function of solar activity 

cost = .3 = amplitude of temperature variation 
cost = 2.2 
cost = 3.0 
| <£+£© |/2 

\4>- Sq \/2 

I H + 0 + p sin (H+7) | 
latitude of the satellite 
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*0 

H 


T 

7 


declination of the sun 
hour angle of the sun 

lag of the temperature maximum with respect to the sun * s 
culmination = cost = -37° 

asymmetry of temperature distribution = cost = +6° 
location of temperature distribution asymmetry = cost = +43° 


Figures 9a, 9b, 9c and 9d show the density and the local exospheric temperature 
variations vs. time (1 orbit) in two different seasonal situations: Spring 
equinox and Summer solstice. 


The deceleration of the i^-mass m A due to drag is given by: 


Fd.i = -jffiPiKIVi 


(2.3.6) 


where : 


A 

Ai 

Cd,i 

Pi 

hi 

f 

Vi 

Ai 

R© 


ballistic coefficient of the i th -mass = — 2 

2mi 

cross section of the i th -inass 
drag coefficient of the i th -mass 
p (hi) = atmospheric density at height hi 
height of the i th -mass = Ri-R$ (l-fsin 2 Ai) 
flattening factor of the Earth 

velocity of the i th -mass relative to the atmosphere 
latitude of the i th -mass 
earth* s equatorial radius 


We assume that the atmosphere rotates rigidly with the Earth and, differently 
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from other authors, we take into account the relative motion between the orbit- 
ing frame and the i th -mass mi. The relative velocity is therefore given by: 


Vi — Ri — x Ri 


(2.3.7) 


where «© is the Earth’s angular velocity while R t is the radius vector from the 
Earth’s center to the i th -mass . We have 


Si = R 0 + ?i 

Ri = r i + ft x Ri (2.3.8) 


where R 0 is the radius vector of the origin of the orbiting reference frame, ri 
is the radius vector from that origin to the i^-mass and ft is the orbital rate. 
Hence 


Vi = (ft -ft®) x (R 0 + ?i) + r± 


(2.3.9) 


The aerodynamic forces (per unit mass) are readily computed by substituting 
equation (2.3.9) into equation (2.3.6). 














LOCAL EXOSPHERIC TEMPERATURE 
INC - 28. S deg 
SUMER SOLSTICE 















Page 46 


2.3.2 J 2 Model - 

For satellites other than the geosynchronous ones the zonal harmonics of 
the Earth's gravity field are the major source of orbital perturbations. The 
expression of the potential of the Earth's gravity field for the i^-mass of the 
system, limited to the zonal harmonic terms, is given by: 


v ‘- -i 

f 1 -t (! 

Jn Pno (sin ^i) 



M/ 


where R© is the earth equatorial radius, /i the earth gravitational constant, J n 
are the zonal harmonic coefficients, P no are Legendre polinomials of argument 
sin^i , 8 ± is the latitude of the i^-mass and Ri is the radius vector of the same 
mass. If we limit the expansion of equation (2.3.10) to the J 2 term (J x = 0 in 
the geocentric reference frame) we obtain the expression of the J 2 perturbation 
potential as follows: 


vf 2 


/ \ 3 

_ jx_ /R$\ Jj 

R t \Ri.) 2 


(1 — 3 sin a 5i) 


(2.3.11) 


The gravity perturbation force due to J 2 is therefore given by: 


ff 2 = -grad (vf 2 ) 


(2.3.12) 


Since in functional form vf 2 = vf 3 [R* , sintfj (Ri) ] we have: 


F J2 = 

XiJ 


c>V, 2 

aRi 


avf 2 a ( sin ^i) 
d(sin 8 ± ) dRi 


dRi 

dx 


ij 


j = 1,2,3 
i = 1, . . . , n 


(2.3.13) 


where x u = X L , x i2 = Y± , x i3 = Z L are the components of the radius vector R t of 
the i th -mass in the inertial reference frame. By substituting equation (2.3.11) 
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into equation (2.3.12) we finally get: 

= - f M J 2 4 (1 — 5sin 2 5i) ^ 

F il = ~ I M (1 - 5sin J 5i) ^ (2.3.14) 

F £ = " f ** Jj *• (3-5sin^ ± ) 

Figure 10 shows the components in inertial reference frame of the J 2 gravity 

force over one orbit with an inclination of 28°.5 and an orbital altitude of 450 

km. The components of the J 2 gravity force in inertial reference frame, given 
by equations (2.3.14), must be transformed into orbiting reference frame compo- 
nents in order to be introduced into our dynamics simulation computer code. 
This transformation is performed according to: 

W’}. = [T] (2.3.15) 

where the subscript (o) stands for orbiting reference fame and the subscript (I) 
stands for inertial reference frame. The transformation matrix [T] is given by: 

(2.3.16) 

where I , J and K are the unit vectors of the geocentric inertial reference frame 

while i, j and k are the unit vectors of the orbiting reference frame given by: 
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j = k x i 


(2.3.17) 


where R 0 is the radius vector of the origin of the orbiting reference frame and 
V c is the velocity vector. The orbiting reference frame orbits the Earth at 
constant angular rate 0 = V 0 /R 0 . Eigure 11 shows the components of the J 2 grav- 
ity force in orbiting reference frame for the same case shown in the previous 
figure . 


2.3.3 Thermal Model - 

The thermal environment produces a periodical variation of the tether 
length (contractions and expansions) with consequent variations of the accelera- 
tion levels. The thermal inputs into the tether are: 

A) Solar radiation 

B) Earth albedo 

C) Earth infrared radiation 

D) Aerodynamic heating 

E) Internal heating (i.e., Joule effect...) 


The only cooling process is the radiation emitted by the wire. We have devel- 
oped a new routine, neglecting the thermal inputs D and E, which calculates the 
temperature gradient of the tether per unit length as follows: 


T 


0* ± Q B ± Q* 


|C - 


p cm 


27rroreT 4 


(2.3.18) 






TIME (SEC) 





where : 


r 

<7 

€ 

P 

C 

m 


tether radius 

Boltzman constant 

tether emissivity 

volume density of the tether 

heat capacity of the tether 

tether mass 


For the expression of thermal inputs we have followed the line adopted by SKY- 
HOOK [4] . We included the earth albedo and we derived a more sophisticated 
model for the partial and the total eclipse of the wire [5] . The expressions of 
the Q 1 are as follows: 


Q* 

Q B 

Q c 


2 r a i sun, inc 
a • I SUN • a • r f • cos ( 7 ) 
2 7T r a a IR f T© 


(2.3.19) 


where : 

OC 

a 

A 

f 

7 

<*IR 

T® 


tether abs orbit ivity 

earth albedo (annual average) 

tether area 

^ ^ 0 - sin (2 0)/2 J . * Re 

view factor = — and smp = 

7T R 

sun zenith angle 

infrared tether absorbitivity 

earth temperature 


In the computation of Q B and Q c we neglected the solar declination angle, the 
differences between continental and ocean areas and the dependence with the 
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latitude and seasonal variations. We also assumed that the reflectance of the 
earth is isotropic and diffuse and that the spectral distribution of reflected 
radiation is equivalent to the spectral distribution of the incident radiation. 
We compared the thermal behavior of two different materials, under the same 
environmental conditions: stainless steel and kevlar. Eigure 12 is a sketch of 

the geometry relevant to the thermal effects over the orbit. In Figure 13 and 
14 the tether temperature is plotted vs. time (2 orbits) for a 2 mm diameter 
stainless steel and kevlar tether respectively. The initial temperature is the 
same for the two cases . 

The thermal expansion (or contraction) of a tether segment of length l 0 is 
given by: 

A£ = a T £ 0 AT ma * (2.3.20) 


where a? is the tether thermal expansion coefficient. If we assume a tether 
length of 10 km the expansion/contraction of the stainless steel wire and the 
kevlar wire over one orbit at 450 km altitude are respectively as follows: 


J Af.teel | 


~ 7m 


| A£k.vi»r | ^ 2m 


(2.3.21) 


where we have assumed a Tiat ®oi = 1.2xl0‘ 5 0 K _1 and a T>kev i ar = -2xl0' 6 °K‘ 1 while the 
maximum AT have been derived from the previous runs as follows: | AT steo i* X | ~ 58°K 
and | AT kev i,r X | ^ 


95°K . 



Figure 12. 
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2.4 Concluding Remarks 

A survey of various out-of -plane libration control laws has been carried 
out. Consequently a yo-yo control law with amplitude of the tether length vari- 
ation proportional to the amplitude of the out-of -plane libration has been se- 
lected. This control law provides good damping when applied to a (theoretical) 
two-dimensional system. In the actual 3-dimensional 4-mass tethered system, 
however, energy is transferred to the least damped degrees of freedom (the out- 
of-plane lateral deflections are still undamped in the present simulations) in 
such a way as to decrease the effectiveness of the algorithm for out-of-plane 
libration control. The addition of damping algorithms for the out-of-plane lat- 
eral deflections is therefore necessary (these damping algorithms will be devel- 
oped in the next reporting period) . 

A completely new package of major environmental disturbances relevant to 
the altitude range of the Space Station has been also developed. This package 
consists of: a new drag model based on Jacchia's 1977 model, a thermal model of 
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the wire and a model of the J 2 gravity perturbation. 


3.0 PROBLEMS ENCOUNTERED DURING REPORTING PERIOD 


None 


4.0 ACTIVITY PLANNED FOR THE NEXT REPORTING PERIOD 


In the next reporting period the algorithms for damping the out-of-plane 
lateral oscillations will be developed. The investigation of the forced dynam- 
ics of the 4-mass tethered system during the station-keeping phase will also be 


carried out. 



